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As  part  of  our  ongoing  research  into  the  clear  daytime  sky’s  visible  structure,  we  analyze  over  1500  skylight 
spectra  measured  during  a  seven-month  period  in  Granada,  Spain.  We  use  spectral  radiances  measured 
within  3°  fields  of  view  (FOV’s)  to  define  colorimetric  characteristics  along  four  sky  meridians:  the  solar  me¬ 
ridian  and  three  meridians  at  azimuths  of  45°,  90°,  and  315°  relative  to  it.  The  resulting  clear-sky  chroma- 
ticities  in  44  different  view  directions  (1)  are  close  to  but  do  not  coincide  with  the  CIE  daylight  locus,  (2)  form 
V-shaped  meridional  chromaticity  curves  along  it  (as  expected  from  theory),  and  (3)  have  correlated  color  tem¬ 
peratures  (CCT’s)  ranging  from  3800  Kto^K.  We  also  routinely  observe  that  sky  color  and  luminance  are 
asymmetric  about  the  solar  meridian,  usually  perceptibly  so.  A  principal-components  analysis  shows  that 
three  vectors  are  required  for  accurate  clear-sky  colorimetry,  whereas  six  are  needed  for  spectral  analyses. 
©  2001  Optical  Society  of  America 
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1.  INTRODUCTION 

Although  many  researchers  have  studied  the  angular  dis¬ 
tribution  of  visible-wavelength  radiance  across  the  clear 
daytime  sky,  few  have  analyzed  its  color  variations.  A 
1985  article  by  Bohren  and  Fraser1  corrected  some  com¬ 
mon  misconceptions  about  the  theory  of  sky  colors  (for  ex¬ 
ample,  Bohren  and  Fraser  note  that  uncritical  use  of  Ray¬ 
leigh’s  scattering  theory  leads  to  the  false  conundrum 
that  the  clear  daytime  sky  should  be  violet),  but  little  ex¬ 
perimental  work  has  followed  their  lead.  In  fact,  nearly 
all  subsequent  research  has  been  devoted  to  developing 
empirical  and  theoretical  models  of  clear-sky  radiance 
and  luminance  across  the  sky  dome. 2-6  By  contrast,  rela¬ 
tively  few  studies  have  dealt  with  the  angular  distribu¬ 
tion  of  clear-sky  spectra  and  chromaticities.7  10 

Yet  skylight’s  spectral  and  color  distributions  are  often 
as  important  in  theoretical  and  applied  work  as  are  its  ra¬ 
diance  patterns.  For  example,  skylight  spectra  can  dic¬ 
tate  the  design  and  performance  of  spectrally  selective 
solar-energy  devices  such  as  photovoltaic  systems.  Im¬ 
proving  our  knowledge  of  clear-sky  spectra  is  valuable  in 
many  research  areas,  ranging  from  aircraft  flight  simula¬ 
tion  and  daylighting  predictions  to  developing  radiation 
models  for  predicting  climate  change  and  its  biological  ef¬ 
fects.  Here  we  present  one  contribution  to  the  larger 
goal  of  developing  a  modern,  truly  worldwide  basis  for 
measuring  and  analyzing  the  clear  sky’s  colors. 

What  is  the  current  state  of  clear-sky  color  modeling 
and  measurement?  In  1991  Sekine  modeled  skylight  col¬ 
ors,  considering  separately  the  effects  of  aerosol  optical 

0740-3232/2001/020412-09$15.00 


depth,  solar  elevation,  reflected  surface  light,  and  scatter¬ 
ing  angle.11  Several  years  later,  one  of  us  used  digital 
image  analysis12  to  measure  daytime  and  twilight  clear- 
sky  chromaticities  along  meridians  near  the  horizon.13 
Among  other  points,  this  study  showed  that  near-horizon 
chromaticity  curves  need  not  parallel  the  Planckian  (or 
blackbody)  locus.  In  1999,  Chain  et  al. 10  found  an  em¬ 
pirical  relationship  between  correlated  color  temperature 
(CCT)  and  luminance  under  clear  skies.  This  develop¬ 
ment  let  them  use  the  CIE  method14  to  derive  skylight 
spectra  from  luminances  alone. 

Skylight  spectra  have  been  measured  within  narrow 
fields  of  view  (FOV’s)  in  past  decades,15"28  but  the  instru¬ 
ments  available  then  ruled  out  quick,  systematic  mapping 
of  the  entire  sky.  We  now  have  the  luxury  of  fast, 
narrow-FOV  spectroradiometers  that  let  us  compile  such 
quasi-instantaneous  skylight  maps.  Another  reason  to 
revisit  sky  color  is  that  older  papers’  terminology  is  itself 
confusing,  sometimes  not  distinguishing  between  “day¬ 
light”  and  “skylight.”  To  avoid  this  ambiguity,  we  follow 
convention  and  distinguish  between  (1)  hemispheric  day¬ 
light  that  can  include  direct  sunlight  and  (2)  narrow-FOV 
skylight  that  always  excludes  direct  sunlight.  In  fact, 
chromaticity  curves  for  daylight  irradiances  and  skylight 
radiances  usually  will  differ  at  a  given  place  (consider, 
say,  daylight  and  skylight  chromaticity  curves  measured 
for  one  hour  at  fixed  detector  orientations).  This  is  true 
because  skylight,  unlike  daylight,  is  measured  within 
small  FOV’s  and  in  many  different  directions.  Other  pa¬ 
rameters  that  affect  sky  color  are  the  observing  site’s  lo- 
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cation  and  altitude,  the  sky  state,  and  the  length  of  time 
required  for  each  measurement.  We  discuss  these  param¬ 
eters  for  our  site  in  the  next  section. 

In  Sections  3  and  4  we  analyze  the  angular  and  short¬ 
term  temporal  trends  in  skylight  colors  calculated  from 
our  measured  narrow-FOV  spectral  radiances.  We  ac¬ 
quired  these  spectra  during  a  seven-month  period  at  one 
site  for  a  wide  variety  of  solar  elevations  h  0  and  view  di¬ 
rections.  The  resulting  chromaticities  are  considered  not 
only  en  masse,  but  also  as  individual  points  that  form 
chromaticity  curves  as  functions  of  view-elevation  angle  h 
along  particular  sky  meridians  (h  is  positive  above  the  as¬ 
tronomical  horizon).  In  addition,  we  use  principal  value 
decompositions  to  examine  the  correlations  among  our 
skylight  spectra.  These  correlations  tell  us  how  many 
eigenvectors  are  required  in  linear  skylight  models  to 
make  our  color  identification  and  recognition  reliable.29,30 


2.  EXPERIMENTAL  PROCEDURE 

To  measure  spectral  radiances  across  the  sky  dome,  we 
used  a  LI-COR  model  LI- 1800  spectroradiometer  that  we 
coupled  via  fiber-optic  cable  to  a  small  refracting  tele¬ 
scope  with  a  3°  (0.00215  sr)  FOV.31  The  LI-1800  radiom¬ 
eter’s  optical  characteristics  have  been  described  in  detail 
by  other  researchers.32-34  In  our  work  the  LI-1800’s 
spectral  range  was  limited  to  the  visible  (380-780  nm), 
and  its  wavelength  resolution  was  5  nm.  Because  the  LI- 
1800  takes  12  s  to  acquire  a  visible  spectrum,  practically 
speaking  we  could  measure  as  many  as  3  spectra/min. 
We  fixed  the  viewing  telescope  to  a  pyrheliometer  mount 
that  let  us  aim  the  telescope  manually  in  the  desired  view 
direction.  Uncertainties  in  setting  h  and  relative  azi¬ 
muth  4>re\  were  0.1°  and  1°,  respectively. 

Figure  1  shows  the  h,  4>rc\  at  which  we  measured  sky¬ 
light  spectra.  We  chose  four  meridians:  the  solar  me¬ 
ridian  (or  the  clear-sky  principal  plane)  with  </>rel  =  0°, 


relative  azimuth  (b  . 

0°  rel 


Fig.  1.  Azimuths  relative  to  the  sun  4>rA  and  view  elevations  h 
above  the  astronomical  horizon  at  which  we  measured  Granada 
skylight  spectra. 


and  three  meridians  at  <b,v\  of  45°,  90°,  and  315°  angles  to 
it.  Following  convention,  we  make  d>rei  increase  clock¬ 
wise  from  the  sun’s  azimuth  (i.e.,  fArcl  increases  in  the 
same  sense  as  compass  direction).  To  simplify  our  nota¬ 
tion,  we  denote  view  elevations  beyond  the  zenith  as  hav¬ 
ing  the  original  4>le\  and  h  >  90° .  Thus  in  our  notation 
=  45°,  h  =  105°  is  equivalent  to  d>rei  =  225°,  h 
=  75°  (see  Fig.  1).  Along  each  sky  meridian,  we  made 
measurements  in  15°  steps  for  15°  =£  h  =s  165°.  As  a  re¬ 
sult,  each  set  of  all-sky  measurements  consisted  of  44  ob¬ 
servations  that  required  20  min  to  complete.  At  the  be¬ 
ginning  and  middle  of  each  set,  we  calculated  the  sun’s 
position  and  calibrated  the  telescope’s  (/jro[  by  centering 
the  sun  in  its  FOV. 

Our  measurement  site  is  in  Granada,  Spain  (latitude 
37°  ll'N,  longitude  3°  35'  W,  altitude  680  m),  a  city  that 
has  no  significant  sources  of  industrial  aerosols.  Our  vi¬ 
sual  prerequisite  for  taking  measurements  was  that  the 
sky  must  appear  to  be  cloudless.  Between  February  and 
August  1998  we  measured  1567  skylight  spectra  that  met 
this  requirement,  and  their  associated  h0  ranged  from 
-0.6°  to  62.2°. 

3.  COLORIMETRIC  CHARACTERISTICS 

A.  Chromaticity  Coordinates 

Our  spectra  include  view  elevations  and  solar  elevations 
that  yield  chromaticities  typical  of  the  blues  observed  in 
clear  daytime  skies.  Thus  none  of  our  measurements  re¬ 
flect  the  far  wider  chromaticity  gamut  seen  in  the  twilight 
sky,13  a  phenomenon  that  we  are  currently  investigating. 

Figure  2  shows  the  CIE  1931  chromaticities  of  all 
Granada  clear  skylight  measurements  overlaid  with  both 
the  Planckian  locus  and  the  CIE  daylight  locus.35  Figure 
2  also  includes  our  Granada  skylight  locus,  which  is  above 
and  nearly  parallel  to  the  CIE  locus  (i.e.,  toward  the 
greens).  A  least-squares  fit  to  our  observed  chromatici¬ 
ties  gives  the  Granada  skylight  locus  as 

y  =  -0.24770  +  2.72203x  -  2.77935x2.  (1) 

The  chromaticity  gamut  of  our  experimental  clear-sky 
measurements  is  broader  than  earlier  ones,15-28  despite 
the  shorter  time  period  of  our  study.  None  of  our  sky¬ 
light  spectra  yield  chromaticities  below  the  Planckian  lo¬ 
cus,  and  most  chromaticities  (  ~  90%)  are  above  the  CIE 
locus.  For  CCT’s  <  7000  K,  the  Granada  skylight  locus 
is  closer  to  the  greens  than  any  comparable  daylight  and 
skylight  chromaticities  published  thus  far.  Even  though 
Granada’s  skylight  chromaticities  have  a  large  CCT 
range,  our  locus,  Eq.  (1),  describes  them  quite  well.  In 
fact,  the  square  of  the  correlation  coefficient  for  Eq.  (1)  is 
R 2  =  0.967,  meaning  that  our  chromaticities  have  negli¬ 
gible  scatter  in  the  green-pink  directions.  In  fact,  the 
95%  confidence  interval  in  chromaticity  distance  about 
the  locus  of  Eq.  (1)  ranges  from  a  maximum  of  0.001  for 
CCT  >  4000  K  to  a  maximum  of  0.005  for  CCT  <  4000  K 
(because  there  are  fewer  measurements  at  these  smaller 
CCT’s). 

Would  our  skylight  chromaticities  exhibit  more  scatter 
if  we  measure  them  within  a  much  smaller  FOV  (say, 
0.25°)?  To  answer  this  question,  consider  its  corollary: 
Do  Granada  daylight  chromaticities  from  hemispheric 
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Fig.  2.  CIE  1931  chromaticities  of  our  1567  Granada  clear  skylight  spectra  (open  circles),  overlaid  with  the  CIE  daylight  locus 
(dashed  curve),  the  Planckian  locus  (curve  with  open  squares),  and  the  Granada  clear  skylight  locus  [curve  with  open  triangles;  Eq.  (1)]. 
The  inset  shows  the  entire  CIE  1931  diagram  and  Planckian  locus. 


skylight  irradiances  (i.e.,  daylight  minus  direct  sunlight) 
show  significantly  less  scatter  than  our  3°-FOV  skylight 
chromaticities?  When  we  compare  the  two  kinds  of  clear- 
sky  chromaticities,  in  fact  we  find  no  statistically  signifi¬ 
cant  decrease  in  the  scatter  of  daylight  chromaticities, 
even  though  they  are  derived  from  a  cosine-corrected 
hemispheric  FOV.  This  suggests  that  if  we  measured 
skylight  spectra  by  using  a  much  smaller  FOV,  the  scat¬ 
ter  in  our  results  would  not  change  in  any  meaningful 
way.  Of  course,  this  is  very  different  from  our  statement 
above  that  daylight  and  skylight  chromaticity  curves  gen¬ 
erally  will  differ  at  a  given  place. 

B.  Correlated  Color  Temperatures 

Numerous  researchers  have  shown  that  CCT  provides  ac¬ 
curate  estimates  of  daylight’s  visible-wavelength  power 
spectra,15-28  and  this  makes  CCT  a  useful  shorthand  for 
specifying  the  colorimetric  and  spectral  characteristics  of 
daylight.  The  CIE  exploits  this  relationship  by  using 


inverse  CCT  (MK'1) 

Fig.  3.  Histogram  of  inverse  CCT  (in  reciprocal  megaKelvins) 
corresponding  to  our  Granada  skylight  chromaticities.  Some 
low-frequency  outliers  are  excluded. 


CCT’s  from  4000  K  to  25,000  K  to  estimate  the  relative 
spectral  irradiances  Ek  of  typical  daylight  phases  at  the 
Earth’s  surface.36 

Our  spectroradiometer’s  software  is  not  designed  to  cal¬ 
culate  CCT’s  >  105  K,  chiefly  because  it  relies  on  older 
techniques  for  approximating  CCT.  This  design  limitation 
sometimes  posed  problems  for  us  in  measuring  skylight 
CCT’s,  since  these  sometimes  exceed  105  K.  In  a  recent 
paper  we  described  a  simple,  accurate  CCT  algorithm 
that  avoids  this  limitation.37  For  the  Granada  skylight 
chromaticities,  our  algorithm  calculated  CCT’s  ranging 
from  3804  K  to  more  than  105  K  (six  chromaticities  had 
CCT’s  >  105  K).  Among  our  1567  skylight  chromatici¬ 
ties,  only  two  were  too  far  from  the  Planckian  locus  to 
yield  spectrally  meaningful  CCT’s. 

When  analyzing  daylight  and  skylight,  often  it  is  more 
convenient  to  use  inverse  CCT  rather  than  CCT  itself. 
The  unit  of  inverse  CCT  is  the  reciprocal  megaKelvin 
(106/CCT),  denoted  by  the  symbol  MK  1  (historically 
known  as  the  mired).38  Reciprocal  megaKelvins  produce 
a  uniform  scale  that  better  describes  human  color  sensi¬ 
tivity  than  does  CCT  proper.  On  this  inverse  scale,  the 
CCT  range  of  our  Granada  clear-skylight  data  is  1.7- 
263.8  MK-1. 

When  this  data  is  plotted  on  a  histogram  (Fig.  3),  the 
most  frequent  inverse  CCT  interval  is  70-75  MK-1.  Fig¬ 
ure  3  can  be  compared  only  imperfectly  with  earlier  work, 
primarily  because  authors  present  their  results  in  widely 
differing  ways  and  because  their  analyses  sometimes  do 
not  distinguish  between  skylight  and  daylight.15-28  That 
caveat  aside,  however,  Fig.  3  does  differ  substantially 
from  its  skylight-only  predecessors. 15'20,21,25  Perhaps  the 
most  significant  departure  in  Fig.  3  is  its  near-normal  dis¬ 
tribution  of  inverse  CCT,  a  pattern  observed  only  once 
before.25 

C.  Meridional  Chromaticity  Curves 

Figures  4-7  show  the  CIE  1931  chromaticities  that  we 
measured  at  different  h  along  the  solar  meridian  and  the 
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meridian  perpendicular  to  it  (  <bni  =  0°,  90°,  respectively) 
at  four  different  solar  elevations  (h0  ~  60°,  20°,  10°,  1°) 
on  three  different  days.  Because  topography  did  not 
block  our  view  of  the  sky  at  h  —  10°  and  h  =  5°  on  two 
days,  we  also  include  those  skylight  chromaticities  in 
Figs.  5  and  7.  To  aid  comparisons  among  Figs.  4-7,  we 
include  the  entire  1931  CIE  diagram  as  an  inset  in  each. 
Within  each  inset,  we  draw  a  small  box  to  show  the  region 
plotted  in  the  figure.  The  chromaticity  curves  in  Figs. 
4-7  are  typical  of  those  measured  on  many  other  days. 

Not  surprisingly,  as  h0  decreases,  the  chromaticity 
gamut  measured  along  the  solar  meridian  increases,  ex¬ 
panding  toward  the  yellow-orange  near  the  sun  [see  espe¬ 
cially  Fig.  7(a)].  This  trend  is  easy  to  interpret 
physically:  Slant-path  optical  depth  increases  with  de¬ 
creasing  h  and  h0,  and  the  resulting  increase  in  direct- 
beam  scattering  makes  the  solar  aureole  yellower  at  a 
given  angular  radius  (or  scattering  angle  1F). 

However,  for  scans  along  </>rel  =  90°,  the  Granada  sky¬ 
light  chromaticity  gamuts  (1)  depend  very  little  on  h0  and 
(2)  are  consistently  smaller  than  those  along  the  solar  me- 


I  1  I  1  I  •  I 
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Fig.  4.  (a)  Chromaticity  curve  for  Granada  clear  skylight 

along  tf>lA  =0°  (the  clear-sky  principal  plane)  when  solar  eleva¬ 
tion  h0  =  60°  on  August  19,  1998.  Note  that  view  elevations  be¬ 
yond  the  zenith  are  indicated  by  h  >  90°.  The  dotted  line  is  the 
CIE  daylight  locus,  while  the  inset  shows  the  entire  CIE  1931 
diagram  and  Planckian  locus.  Figures  5-7  reuse  Fig.  4’s  nota¬ 
tion.  (b)  Chromaticity  curve  for  clear  skylight  along  </>rel 
=  90°  when  h0  =  60°  (August  19,  1998). 
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Fig.  5.  (a)  Chromaticity  curve  for  clear  skylight  along  </>rel 

=  0°when/i0  =  19°  (August  6,  1998).  (b)  Chromaticity  curve 
for  clear  skylight  along  </>rel  =  90°  when  h0  =  19°  (August  6, 
1998). 

ridian.  Only  for  h0  >  30°  [see  Fig.  4(b)]  does  the  (bre[ 
=  90°  gamut  shrink  appreciably.  In  fact,  the  chroma¬ 
ticities  in  Fig.  4(b)  span  a  Cartesian  distance  of  ~9  just- 
noticeable  differences39  (JND’s)  while  their  solar- 
meridian  counterparts  in  Fig.  4(a)  cover  a  distance  more 
than  twice  as  large.  For  h0  <  30°,  chromaticity  gamuts 
along  the  <b„.\  =  90°  meridian  are  nearly  independent  of 
sun  elevation. 

As  Sekine  shows,  chromaticity  curves  along  the  solar 
meridian  are  V-shaped  (possibly  with  a  zigzag  near  h0 
when  the  sun  is  high  in  the  sky),  and  they  nearly  parallel 
the  Planckian  locus  on  its  green  side.11  The  antisolar 
half  of  these  curves  [i.e.,  <£rel  =  0°,  h  >  90°;  see  Fig.  4(a)] 
is  closer  to  the  greens  than  is  the  solar  half.  The  same 
V-shaped  pattern  is  evident  in  our  measured  chromaticity 
curves  for  </>rel  =  90°  [Figs.  4(b),  5(b),  6(b),  and  7(b)],  a 
relative  azimuth  that  Sekine  did  not  model.  Our  mea¬ 
surements  also  clearly  show  that  the  4>re\  ~  0°  chromatic¬ 
ity  curves  nearly  coincide  with  the  CIE  locus  (dashed  line 
in  Figs.  4-7),  whereas  their  (b„.\  =  90°  counterparts  lie 
some  distance  from  it.  This  departure  from  the  CIE  lo¬ 
cus  is  largest  at  low  CCT’s  (<10,000  K),  where  the  maxi¬ 
mum  chromaticity  distance  just  exceeds  4  JND’s. 

One  of  the  longest-standing  assumptions  about  cloud¬ 
less  skies  is  that  their  chromaticity  and  luminance  distri- 
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butions  are  symmetric  about  the  principal  plane,  includ¬ 
ing  along  the  </irei  =  90°  meridian.  However,  our 
measurements  show  something  quite  different  in  Fig.  8. 
To  generate  it,  we  calculate  chromaticity  distances  be¬ 
tween  pairs  of  skylight  colors  measured  at  the  same  view 
elevation  h  along  d>rel  =  90°  (say,  at  h  —  60°,  120°).  We 
then  normalize  each  chromaticity  distance  by  the  appro¬ 
priate  JND  and  plot  the  ratios  at  the  appropriate  h. 
When  h0  —  60°,  the  maximum  chromaticity  distance 
along  4>re\  ~  90°  is  ~2  JND  at  h  =  60°  (open  squares  in 
Fig.  8;  see  right  ordinate).  At  lower  h0  the  maximum 
chromaticity  distance  exceeds  4  JND’s  at  most  h  (open, 
upright  triangles  in  Fig.  8).  Thus,  regardless  of  sun  el¬ 
evation,  color  asymmetries  at  the  same  h  on  either  side  of 
the  principal  plane  are  often  perceptible.  In  other  words, 
if  we  simultaneously  compared  adjacent  samples  of  two 
such  colors,  we  could  distinguish  between  them. 

Not  unexpectedly,  we  find  similar  luminance  asymme¬ 
tries  about  the  principal  plane.  In  Fig.  8  the  largest  lu¬ 
minance  differences  (~20-25%;  see  left  ordinate)  occur  at 
the  lowest  solar  elevation  (h 0  =  0.4°),  although  these  are 
paired  with  fairly  small  chromaticity  differences.  For 
the  remaining  h0,  luminance  differences  are  <8%. 
Nonetheless,  if  we  take  a  2%  luminance  difference  as  a 
typical  threshold  contrast  for  photopic  vision,  then  nearly 


Fig.  6.  (a)  Chromaticity  curve  for  clear  skylight  along  </>rel 

=  0°  when  h0  =  11°  (August  19,  1998).  (b)  Chromaticity 

curve  for  clear  skylight  along  </>rel  =  90°  when  h0  =  11°  (August 
19,1998). 


x  CIE  1931 


(a) 


Fig.  7.  (a)  Chromaticity  curve  for  clear  skylight  along  tf. >rel 

=  0°  when  h0  =  1.6°  (August  5,  1998).  (b)  Chromaticity 

curve  for  clear  skylight  along  </>rel  =  90°  when  h0  =  1.6°  (August 
5,1998). 


Fig.  8.  Relative  luminance  differences  (solid  symbols)  and 
chromaticity  differences  in  JND  units  (open  symbols)  as  a  func¬ 
tion  of  h  along  </>rel  =  90°  for  the  four  h0  shown  in  Figs.  4(b),  5(b), 
6(b),  and  7(b).  Solar  elevation  symbols  are  squares  for  h0 
=  60°,  circles  for  h0  =  18.5°,  upright  triangles  for  h0  =  10°, 
and  inverted  triangles  for  h0  =  0.4°. 

all  of  the  clear  sky’s  luminance  asymmetries  about  the 
principal  plane  are  perceptible  (naturally,  we  rarely  make 
the  required  side-by-side  comparisons  of  skylight 
samples).  Only  for  h0  =  18.5°  are  the  Granada  skylight 
luminance  differences  consistently  subthreshold.  Our 
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measurements  suggest  that  such  skylight  luminance  and 
color  asymmetries  are  the  norm,  rather  than  the  excep¬ 
tion,  in  real  atmospheres. 

Some  of  this  asymmetry  might  be  ascribed  to  the  time 
required  to  make  our  meridional  measurements.  Yet 
during  the  4  min  needed  for  each  scan  h0  changes  by  only 
±0.5°  at  high  solar  elevations  and  ±1°  at  low  ones.  Be¬ 
cause  our  uncertainty  in  measuring  </>rel  is  of  comparable 
magnitude,  we  believe  that  changes  in  solar  position  are 
unlikely  to  cause  the  observed  asymmetries  in  the  lumi¬ 
nance  and  chromaticity  of  Granada’s  clear  skies.  This 
claim  is  supported  by  clear-sky  spectra  that  we  measured 
with  another  fast-scanning  spectroradiometer.  Although 
we  acquired  spectra  at  matching  h  and  <pre\  within  a  few 
seconds  of  each  other,  we  still  found  perceptible  differ¬ 
ences  in  skylight  color  and  luminance  about  the  principal 
plane. 

One  plausible  explanation  for  such  asymmetries  is  that 
aerosol  content  varies  with  tbrc\ .  In  Granada  this  might 
be  caused  by  enhanced  aerosol  loading  from  farmland 
west  of  our  measurement  site.  The  skylight  asymmetries 
are  less  likely  to  be  caused  by  changes  in  surface  reflec¬ 
tance  spectra,  because  these  do  not  vary  in  any  system¬ 
atic  way  within  40  km  of  our  site.  Although  most  extant 
models  of  atmospheric  scattering  assume  symmetry  about 
the  principal  plane,2"6  our  results  show  that  this  assump¬ 
tion  is  only  a  first  approximation,  not  an  inviolable  fact. 
Thus,  to  model  clear-sky  optics  more  realistically,  we 
should  add  the  asymmetries  in  color  and  luminance  re¬ 
ported  here  to  the  asymmetries  in  polarization  that  we 
measured  earlier.40  Although  we  have  not  yet  explicitly 
modeled  such  skylight  asymmetries,  recognizing  their  ex¬ 
istence  is  a  necessary  first  step  in  doing  so. 

We  complete  our  colorimetric  analysis  in  Figs.  9-12  by 
showing  chromaticity  curves  along  the  4>re\  =  45°,  315° 
meridians  for  h0  =  58°,  52°,  42°,  6°.  In  each  figure  the 
meridional  scans  have  similar  shapes,  since  both  are  45° 
from  the  principal  plane.  However,  chromaticities  can 


Fig.  9.  Chromaticity  curves  for  Granada  clear  skylight  along 
<f>Te\  =  45°,  315°  when  solar  elevation  h0  =  58°  on  August  19, 
1998.  Chromaticities  along  </>rel  =  45°  are  denoted  with  filled 
squares;  those  along  </>rel  =  315°,  with  open  circles.  As  in 
Figs.  4-7,  view  elevations  beyond  the  zenith  at  each  </>rel  are  in¬ 
dicated  by  h  >  90° .  The  dotted  line  is  the  CIE  daylight  locus, 
while  the  inset  shows  the  entire  CIE  1931  diagram  and  Planck- 
ian  locus.  Figures  10-12  reuse  Fig.  9’s  notation. 
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4>tA  =  45°,  315°  when  h0  =  52°  (August  20,  1998). 


Fig.  11.  Chromaticity  curves  for  Granada  clear  skylight  along 
4>le\  =  45°,  315°  when  h0  =  42°  (August  25,  1998). 


x  CIE  1931 

Fig.  12.  Chromaticity  curves  for  Granada  clear  skylight  along 
4>re l  =  45°,  315°  when  h0  =  6°  (August  6,  1998). 


differ  markedly  at  the  same  scattering  angle,  such  as  Fig. 
12’s  two  h  =  165°  chromaticities.  At  both  relative  azi¬ 
muths,  the  chromaticity  gamut  increases  with  decreasing 
h0  as  skylight  is  increasingly  reddened  by  direct  sunlight. 
Note,  too,  that  the  positions  of  the  chromaticity  curves 
shift  from  day  to  day,  probably  as  a  result  of  changes  in 
aerosol  optical  depth.11 
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Fig.  13.  Inverse  CCT  as  a  function  of  scattering  angle  'F  for 
all  Granada  skylight  chromaticities.  Solid  squares  mark  the 
mean  inverse  CCT  for  each  5°-wide  bin,  and  error  bars  are  2 
standard  deviations  wide. 

D.  Scattering  Angle 

Our  most  vivid  Granada  skylight  had  a  colorimetric  pu¬ 
rity  of  39.4%,  while  the  average  skylight  purity  was  21.6% 
(measured  with  respect  to  sunlight’s  chromaticity  outside 
the  atmosphere41).  Nearly  33%  of  our  skylight  spectra 
had  purities  >25%,  and  the  modal  dominant  wavelength 
for  all  spectra  was  478-479  nm.  Compared  with  the  the¬ 
oretical  Rayleigh  purity  of  41%  for  a  single-scattering, 
purely  molecular  atmosphere,  our  measured  purities  may 
seem  small  indeed,  yet  they  are  similar  to  those  we  have 
measured  elsewhere.13  As  simple  theory  predicts,  we 
found  that  the  purest  blues  at  Granada  (1)  occur  near  4' 
~  90°  (rather  than  at  the  zenith,  where  optical  depth  is 
smallest)  and  (2)  vary  daily  in  both  dominant  wavelength 
and  purity,  with  the  large  optical  depths  of  hazy  days  pro¬ 
ducing  the  smallest  maximum  purities.1  However,  as 
noted  below,  the  highest  skylight  purities  occur  not  at  4/' 
=  90°,  but  at  somewhat  larger  scattering  angles. 

In  Fig.  13  we  analyze  the  scattering-angle  dependence 
of  Granada’s  skylight  chromaticities  in  the  form  of  in¬ 
verse  CCT.  We  first  sorted  our  measurements  into  36 
scattering-angle  bins  between  0°  and  180°,  with  each  bin 
containing  inverse  CCT’s  from  a  5°  interval  in  4'.  We 
then  calculated  the  mean  inverse  CCT  and  its  standard 
deviation  for  each  bin.  These  standard  deviations  indi¬ 
cate  the  chromaticity  variability  caused  by  changes  in  (1) 
h0  and  (2)  slant-path  molecular  and  aerosol  optical 
depths.  Because  aerosols  scatter  anisotropically,  Fig.  13 
is  not  exactly  symmetric  about  90°,  as  it  would  be  in  a 
purely  molecular  atmosphere.  In  fact,  aerosols’  greatly 
enhanced  forward  scattering  causes  the  lowest  inverse 
CCT’s  (the  purest  blues)  to  occur  not  at  1F  =  90°,  but  in¬ 
stead  from  90°  <  T  <  110°,  close  to  the  minimum  of  the 
aerosol-scattering  phase  function. 


4.  SPECTRAL  ANALYSIS 

Spectral  analysis  of  our  skylight  measurements  is  also  in¬ 
structive,  and  we  do  this  with  a  principal  value  decompo¬ 
sition  of  the  kind  that  Parkinnen  et  al.  describe.42  Let 
Ee(\)  be  an  experimentally  determined  skylight  spec¬ 
trum.  We  can  reconstruct  it  approximately  by  using  p 
eigenvectors  in  the  equation 


p 

ER{\)  =  2  {Ee(X)\ Y;(k))Vri(\),  (2) 

i  =  1 

where  ER (A.)  is  the  reconstructed  spectrum,  Vj( X )  the  ith 
eigenvector,  and  (|)  denotes  the  inner  product.  To  test 
this  reconstruction’s  accuracy,  we  once  again  use  a  good¬ 
ness  of  fit  coefficient43-45  (GFC)  that  is  based  on  the 
Schwartz’s  inequality.  In  our  work,  this  coefficient  is 


GFC 


2  EE(\j)ER(\j) 


2 


IEe(\j)]2 


[Er^j)]2 


1/2  ■ 


(3) 


In  fact,  the  GFC  is  the  multiple  correlation  coefficient 
R,  the  square  root  of  ER(\)’s  spectral  variance  with  re¬ 
spect  to  the  original  Ee(\).  The  GFC  ranges  from  0  to  1, 
with  1  corresponding  to  an  exact  duplicate  of  Ee(X).  We 
found  that  colorimetrically  accurate  E  R(X)  require  a 
GFC  >  0.995;  what  we  call  “good”  spectral  fit  requires  a 
GFC  >  0.999,  and  GFC  >  0.9999  is  necessary  for 
“excellent”  fit.43’44 

We  plot  the  first  five  eigenvectors  V(X)  for  our  Granada 
skylight  dataset  in  Fig.  14,  and  collectively  they  account 
for  99.9%  of  its  measured  spectral  variance.  The  second 
skylight  V(X)  corresponds  to  the  dataset’s  yellow-blue 
variations,  as  earlier  research  has  shown.18,23  Yet  in  a 
departure  from  earlier  results,  our  third  skylight  V(X) 
does  not  account  for  the  dataset’s  purple-green  varia¬ 
tions.  Instead  this  eigenvector  produces  Granada  chro¬ 
maticity  displacements  that  nearly  parallel  the  CIE  day¬ 
light  locus. 

Although  we  use  the  CIE  method  to  recover  skylight 
spectra  from  their  corresponding  1931  chromaticities,36 
we  update  it  by  using  the  Granada  skylight  locus  [Eq.  (1)] 
and  skylight  eigenvectors  (Fig.  14).  The  CIE  suggests 
using  just  the  first  three  TOO,  which  for  our  data  yields  a 
mean  GFC  of  just  0.9900:  Only  53.3%  of  our  ER(X)  are 
of  “acceptable”  quality,  while  “good”  and  “excellent” 
ER(X)  are  limited  to  only  21.4%  and  1.91%,  respectively. 


400  450  500  550  600  650  700  750 


Wavelength  (nm) 

Fig.  14.  Spectral  distribution  of  the  first  five  eigenvectors  for 
the  Granada  skylight  basis  set  from  380-780  nm.  Solid  curve, 
p  =  1  eigenvector  (mean  vector);  dashed  curve,  p  =  2  eigenvec¬ 
tor;  dotted  curve,  p  =  3  eigenvector;  solid  curve  with  open 
squares,  p  =  4  eigenvector;  dashed  curve  with  solid  circles,  p 
=  5  eigenvector. 


Hernandez-Andres  et  al. 


Fig.  15.  Mean  CIELUV  color  difference  A E*uv  as  a  function 
of  mean  GFC.  Numbers  next  to  each  point  indicate  the  number 
of  eigenvectors  p  necessary  in  Eq.  (2)  to  yield  the  given  GFC. 
Error  bars  are  2  standard  deviations  wide. 

This  limited  spectral  accuracy  arises  because  the  CIE 
method  permits  reconstruction  of  exact  ER (\)  only  if  we 
use  chromaticities  exactly  on  the  given  locus  [Eq.  (1), 
here].  Clearly  very  few  of  our  measured  chromaticities 
meet  this  criterion,  and  this  reduces  the  mean  spectral  ac¬ 
curacy  of  our  Er{\). 

What  happens  if  we  step  outside  the  CIE  recommenda¬ 
tions  and  test  different  p  values  in  Eq.  (2)?  Now  we  cal¬ 
culate  the  resulting  GFC  values  and  CIELUV  color  differ¬ 
ences  A.E*„„  between  the  original  and  reconstructed 
spectra.  As  we  see  in  Fig.  15  we  must  use  at  least  the 
first  three  V(\)  in  order  to  get  “good”  mean  quality, 
whereas  “excellent”  ER(\)  require  live  V(\).  For  more 
than  95%  of  our  E  lt(K)  to  be  of  “acceptable,”  “good,”  or 
“excellent”  quality,  we  require  2,  4,  or  6  V(\),  respec¬ 
tively.  Although  in  principle  metamerism  can  decouple 
colorimetric  accuracy  from  spectral  accuracy,  Fig.  15  indi¬ 
cates  that  in  principal  our  mean  A E*uv  depends  strongly 
on  the  average  spectral  GFC. 

Figure  15  also  shows  the  dependence  of  mean  A E*uv  on 
the  number  of  V(X)  used.  As p  increases,  colorimetric  er¬ 
rors  decrease  sharply,  especially  for  the  first  three  eigen¬ 
vectors.  In  technical  or  industrial  applications,  three  or 
five  A E*uv  units  are  often  taken  to  be  a  JND.  Figure  15 
shows  that,  by  using  p  =  3  in  Eq.  (2),  we  ensure  not  only 
a  mean  A E*uv  <  3  units,  but  also  that  its  associated 
standard  deviation  is  <  3  units.  Our  analysis  shows,  in 
agreement  with  CIE  recommendations,  that  in  fact  three 
VOO  are  required  to  reconstruct  a  skylight  spectrum  that 
is  colorimetrically  and  photometrically  indistinguishable 
from  its  original  Ee(X). 


5.  CONCLUSIONS 

Among  the  many  new  observations  stemming  from  our 
skylight  measurement  campaign,  several  stand  out. 
First,  the  CIE  daylight  locus  only  approximates  Grana¬ 
da’s  great  variety  of  clear  skylight,  which  has  CCTs  rang¬ 
ing  from  3804  K  to  °°  K  and  inverse  CCT  peaks  at  70-75 
MK  .  Second,  our  new  Granada  locus  [Eq.  (1)]  more  ac¬ 
curately  measures  local  skylight  variability  than  does  the 
CIE  locus.  In  particular,  our  locus  is  farther  from  the 
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greens  than  any  other  published  daylight  or  skylight  loci. 
This  in  turn  suggests  the  value  of  developing  different 
skylight  loci  for  different  types  of  sites  rather  than  relying 
on  just  one  locus  worldwide. 

Third,  chromaticity  curves  along  the  4>re l  =  0°,  90°  me¬ 
ridians  (the  principal  plane  and  its  perpendicular)  form 
distinct  V  shapes  that  nearly  coincide  with  the  CIE  locus. 
Fourth,  at  all  relative  azimuths  we  found  distinct  asym¬ 
metries  in  color  and  luminance  at  geometrically  equal  h, 
d>rel .  These  differences  usually  are  suprathreshold,  and 
they  can  be  as  large  as  4-5  JND  or  8%  in  luminance.  Con¬ 
trary  to  the  simplifying  assumptions  made  in  many  clear- 
sky  scattering  models,  such  asymmetries  are  routinely 
perceptible  in  the  real  atmosphere.  Fifth,  our  principal 
value  decomposition  shows  that  three  eigenvectors  (in 
agreement  with  CIE  recommendations)  are  needed  to  re¬ 
construct  colorimetrically  accurate  skylight  spectra. 
However,  six  eigenvectors  are  required  for  spectrally  ac¬ 
curate  skylight  spectra. 

Although  skylight  chromaticities  and  spectra  have 
been  measured  precisely  for  several  decades,  our  work 
shows  that  new  measurements  can  still  spur  valuable 
new  theoretical  and  practical  insights.  For  example,  our 
quantification  of  clear-sky  color  and  luminance  asymme¬ 
tries  suggests  how  much  variability  a  truly  comprehen¬ 
sive  visible-wavelength  radiative  transfer  model  must  ac¬ 
commodate.  As  for  practical  applications,  our  work 
easily  lets  us  calculate  (1)  diffuse  sky  spectral  irradiance 
on  sloping  surfaces  and  (2)  skylight  chromaticities  that 
are  angularly  integrated  over  FOV’s  much  larger  than  3°. 
Thus  anytime  or  anywhere  that  skylight  rather  than  day¬ 
light  proper  dominates  the  daytime  sky,  our  new  observa¬ 
tions  can  have  great  utility. 
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